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Abstract 

In dynamic force spectroscopy, a (bio-)molecular complex is subjected to 
a steadily increasing force until the chemical bond breaks. Repeating the 
same experiment many times results in a broad distribution of rupture forces, 
whose quantitative interpretation represents a formidable theoretical chal- 
lenge. In this study we address the situation that more than a single molecu- 
lar bond is involved in one experimental run, giving rise to multiple rupture 
events which are even more difficult to analyze and thus are usually elimi- 
nated as far as possible from the further evaluation of the experimental data. 
We develop and numerically solve a detailed model of a complete dynamic 
force spectroscopy experiment including a possible clustering of molecules 
on the substrate surface, the formation of bonds, their dissociation under 
load, and the post processing of the force extension curves. We show that 
the data, remaining after elimination of obvious multiple rupture events, 
may still contain a considerable number of "hidden" multiple bonds, which 
are experimentally indistinguishable from "true" single bonds, but which 
have considerable effects on the resulting rupture force statistics and its 
consistent theoretical interpretation. 



1 Introduction 

The specific binding of a ligand molecule to a receptor protein is an essential 
functional principle of molecular recognition and information processing in 
biological systems, e.g. in the context of genome replication and transcrip- 
tion, enzymatic activity, initiation of infection and immune response (1). 
Furthermore, their force-induced dissociation is of great interest with respect 
to adhesion and cohesion of any type of biological matter and their selectiv- 
ity can be exploited in various bioanalytical and biomedical devices (j^, y) . 
The invention and continuous refinement of atomic force microscopy (AFM) 
and other related techniques p, yj) has led to a tremendous progress in our 
abilities to directly investigate such processes on the molecular level (jj, M) ■ 
Accordingly, dynamic force spectroscopy has become a standard tool in ex- 
perimental Biophysics to explore and quantify molecular recognition and 
binding properties of (bio-) molecular complexes, such as reaction off-rates, 
binding energy landscapes, affinity ranking of single point mutations etc. 
([l|43|). Typical examples include ligand-receptor compounds like antibody- 
antigen (|6|), protein-DNA (|7|), or supramolecular complexes ([g). In all these 
experiments, one of the participating molecules is connected to a tiny force 
probe and the other to a substrate surface (or an immobilized micron-sized 
bead). By controlling the relative distance s between these objects with 
nanometer precision, an experimentally controlled, steadily increasing force 
F(s) in the Piconewton range can be applied to the compound until the 
chemical bonds between the constituting molecules rupture, see Fig. [T]A. 

Repeating the same experiment many times results in a broad distribu- 
tion of rupture forces, whose interpretation and quantitative analysis rep- 
resents a formidable theoretical challenge, of foremost importance for our 
understanding of the basic biophysical principles at work and for the prac- 
tical evaluation of experimental rupture force data. 

In other words, given the distribution of rupture forces for various pulling 
velocities, a theoretical framework is needed in order to achieve a quantita- 
tive characterization of the probed chemical bond in terms of binding affini- 
ties, force free dissociation rates, and molecular binding energy landscape 
parameters. The predominant general framework in this context is due to 
Evans and Ritchie (l_9i) , describing the rupture of a molecular complex as the 
dissociation of a single chemical bond under the externally applied force in 
terms of a thermally activated escape process with a force dependent decay 
rate. While Evans and Ritchie adopted the phenomenological approach by 
Bell did); important subse quen t refinements concern e.g. the specific force 
dependence of the rate (Ill4l4l). multiple dissociation pathways (jlSl). linker 



length distributions (jld). and contact times (jlTI ). 

Experimentally, in a given force extension curve F{s) it is often quite 
obvious that more than one bond has actually been probed (Figs. [IJ3,C). 
Since such multiple rupture events cannot be quantitatively evaluated within 
the above theories, they are usually eliminated as far as possible from the 
experimental data sets. The key issue of our present work is to take into 
account the possibility that even after eliminating multiple bonds along this 
line, the remaining apparent single bonds may in fact still contain a signif- 
icant number of "hidden" multiple bonds. As we will show, this number 
is indeed far from negligible under usual experimental conditions. While a 
recent experimental work has demonstrated that this might indeed be the 
case (jla). a thorough theoretical study of this point has still been missing so 
far. Such a violation of the basic assumption of the above standard approach 
in this context is not only of conceptual interest, but also may have con- 
siderable practical implications regarding the numerous experimental works 
employing this "standard" theoretical framework for evaluating their data. 

Specifically, within this standard framework the main focus is usually 
put on the dependence of the most probable rupture force /* on the pulling 
velocity v. While the corresponding predictions indeed agree very well with 
most experiments, a more careful analysis of the complete rupture force 
distribution reveals that most experimental data are actually incompatible 



with any of those models ([19|, 1201 ) . in particular the typical "long tails" 
( 171 . Iisl ) of the experimental rupture force distributions, as exemplified by 
Figs. [lfc,F. 

A first, purely phenomenological explanation of those incompatibilities 
in terms of bond- heterogeneities was proposed in (J2Cl ) , amounting to a para- 
metric randomization of the force dependent single bond dissociation rates. 
Our present work provides for the first time a more detailed microscopic 
foundation of this ansatz, since a mixture of "true" single bonds and undis- 
covered multiple bonds is similar to a statistical ensemble of formal single 
bonds with randomized parameters in the force dependent rate (see dis- 
cussion in Sect. 14. 4p . While hidden multiple bonds thus seem to play a 
prominent role in this context, a more detailed analysis, which goes beyond 
the scope of our present paper, however indicates, that there may well ex- 
ist additional microscopic mechanisms which also contribute notably to the 
above mentioned incompatibilities and their explanation in terms of bond- 
heterogeneities. 

Yet another main conclusion of our present paper is that even though 
a considerable fraction of rupture forces are due to (unresolved) multiple 
bonds, the rate parameters can still be extracted from the dependence of 



the most probable rupture force on the pulhng velocity according to the 
above mentioned single bond theories. In other words, the hidden multiple 
bonds heavily contribute to the long tails in Figs. [I](E) and (F), but hardly 
to the main peaks. 

We should emphasize that besides the investigations of ligand-receptor 
dissociation, the general realm of so-called "dynamic force spectroscopy" 
also includes the exploration of protein unfolding, unzipping of DNA strands 



and DNA hairpins, forced rupture of cell adhesion etc., see e.g. (jld . I2ll424l) 
and further references therein. These quite different experimental set ups 
often admit a theoretical treatment quite similar to those mentioned above. 
However, they usually do noi exhibit any kind of incompatibilities e.g. in the 
form of long tails, see e.g. (Ill.ll3l.l23l428l). Our present approach specifically 



and exclusively concerns the case of "dynamic rupture force spectroscopy" 
by AFM, micropipette-based force probes, laser tweezers etc. yHsl). With 
the exception of (1291 ). the problem of the long tails does arise for all those 
rupture force experiments we know of. 

As far as multiple bond models are concerned, we build on several re- 
lated previous studies: A common starting point in many of those works. 






see e.g. (jlOl . |21| . |22| . 13014381) . is the assumption that the applied force is 
equally distributed among all bonds. Closer inspection however shows that 
the distribution of the rupture forces, expected from these models, in general 
exhibits, unlike the experimentally observed ones (Figs. [TJi],F), several well 
separated peaks, which are attributed to single, double ,triple, ... bonds 



(J34l . l38l ). see also Sect. 14.31 below for more details. In other words, the 
assumption of equally distributed forces must be given up. Such a refined 
theoretical treatment has been developed for the first time by Akhremitchev 
and coworkers (|38l ). but their approach still exhibits several other strong 
oversimplifications: First, the number of receptors still cannot be more than 
one, while the number of ligands is at most two. Second, in case of two coex- 
isting bonds, both of them are assumed to rupture simultaneously. Another 
important series of works considering forces which are unequally distributed 



among several coexisting bonds is due to Gao and coworkers ((391, l4d ). For- 
mally, they are quite similar to our present study but physically they focus 
on cell adhesion and do not admit any conclusions with respect to our present 
case. 

Taking into account the very fact that force extension curves with clearly 
visible multiple bond signatures (Figs. [IB,C) are discarded before further 
evaluating those data, is yet another essential difference between our present 
work and all the above mentioned previous multiple bond studies. 

Conceptually, the basic idea of our paper is to establish and analyze a 



simplified but still reasonably realistic model of the entire experimental pro- 
cedure: preparation (functionalization) of the AFM tip and of the sample 
surface; formation of bonds; retraction of the tip from the surface, leading to 
bond dissociation (rupture events); processing and evaluation of the so ob- 
tained data exactly as in the case of real experimental data. The advantage 
of such a "computer experiment" is that the interpretation of the data and 
the final conclusions can be immediately compared with what was actually 
going on in the considered system. 

2 Model 

2.1 Functionalization of tip and sample 

A typical experimental setup we have in mind is exemplified by Fig. [2]A. 
As illustrated by Fig. [2jB, the AFM tip is modeled as a half sphere of 
experimentally realistic radius R = 15 nm, to which ligands are connected 
via linkers of experimentally realistic length L = 30 nm. The number of 
linkers is thus equal to the number of ligands and is denoted by A'lin. The 
vector connecting the tip apex with the immobilization point of linker i 
(i = 1 . . . A'^iin) is denoted as r^^ ™ and the Nn^ immobilization points them- 
selves are randomly sampled according to a uniform distribution on the half 
sphere. In a real experiment, the number A'^iin of linkers attached to the tip 
is itself a random variable with a distribution which depends in a highly 
complicated way on the (chemical) functionalization procedure. Within our 
model, we assume the linker number Nn^ to be given and we will show that 
our conclusions concerning, e.g., the distribution of rupture forces, depend 
very little on the exact value of A'^iin, at least as long as Nn^ > 5 (see Sect. 
[S]). In view of typical grafting densities of, e.g., the PEG linkers used in (jsl), 
a realistic choice which we often adopt in our numerical calculations below 
is Nnn = 10. 

The receptors are randomly immobilized on the sample surface with no 
or negligibly short linkers. Beside the case that they are independently and 
uniformly distributed according to some preset density prec (Sec. [3]), we will 
also discuss the case that the receptors are distributed in clusters of varying 
size (Sec. H]). 

2.2 Force distance cycle 

Usually one refers to a single repetition/run of a pulling experiment as a 
force distance cycle: In a first step, tip and sample are brought into contact 



(at the origin of our coordinate system) for a certain "dwell-time" idweii 
(typically 0.1 — Is), during which receptors and ligands may form bonds. 
Then, the AFM tip is retracted along the z-axis, rapex(0 = (0' 0; •^(t)) (cf. 
Fig. [2]), resulting in a steadily increasing force on the single bonds until 
they rupture. In the force extension curve, bond ruptures induce more or 
less pronounced force dips (Figs. [TJf\-D). As in the real experiment, if more 
than one such dip is clearly visible (Figs. [1|3,C), the corresponding force 
extension curve is discarded. 

Before discussing these steps in more detail in the following subsections, 
we remark that while the AFM tip and thus the linker immobilization points 
r^- remain the same for all cycles (see Fig. [2p), a new random sampling 
of the receptor positions will be performed for each cycle, modeling the fact 
that in a real experiment the sample surface is usually "scanned" either due 
to thermal drift or controlled lateral displacements. 

2.3 Formation of bonds 

We consider N^^ ligands, attached by linkers to the AFM tip, and a surface 
with randomly distributed receptors, having the opportunity to form bonds 
during the dwell time tdweii- The resulting number of bonds thus depends 
on the dynamics of a nanometer-sized object tethered to a polymer chain, 
which itself may interact with other chains, and which is subjected to com- 
plicated boundary conditions (half-sphere on a flat surface) and far from 
equilibrium initial conditions. Instead of theoretically addressing this very 
difficult problem, we invoke the experimental findings in a related situation, 
namely the association of end-grafted chains to uniformly reactive surfaces 



(|41l 443l). As a result of these works, the probability that a bond between 
a tethered molecule and one of its potential binding partners at a distance 
d is actually formed is approximately constant for d < dmax and negligible 
for d > (ijnax; where dmax ~ 10-15 nm for PEG linkers with L ~ 30 nm and 
idwcU ~ Is. Moreover, the number of remaining "unpaired potential binding 
partners" closer than dmax is relatively small. 

We therefore proceed as follows: First, all ligand-receptor pairs with 
a distance smaller than (imax are identified. Then, one of these pairs is 
randoinly chosen to form a bond. For the remaining receptors and ligands 
these two steps are repeated, until no pair with a distance smaller than dmax 
is left. The total number of the so obtained bonds is henceforth denoted as 
N. In our numerical calculation below, our standard choice will be dmax = 
12 nm, while a more detailed exploration of other choices will be discussed 
in Sect. Ol 



We have also tried out other strategies for the bond formation, for exam- 
ple that a ligand preferentially forms a bond with the closest receptor. This, 
however, yielded very similar results to the above mentioned procedure. 

2.4 Force extension curves 

Denoting by r^ " and r^'^'^'^'^ the linker and receptor immobilization points 
of the i-th bond (i = l...A^), the linker's end-to-end vector is (see Fig. [2P) 

, (lin) (rec) , .^ /^x 

Ti = ze^ + r) ' - r) ' = l^n , (1) 

where e^ is the unit vector in z-direction, f j is the unit vector in rj-direction, 
and li is the length (modulus) of rj(z). 

Due to entropic effects, stretching the linker to a length /j requires a 
force of modulus fi which acts in the direction of the end-to-end vector, i.e. 

fi = fiTi . (2) 

A particularly simple model for this force-dependence of the linker length is 
the freely jointed chain (FJC) model, according to which 



hih) = L 



th f ^^^A - ^^'^ 
V ^bT J Xxfi 



(3) 



where Xk is the Kuhn length and ksT the thermal energy. As a typical 
value, obtained by stretching experiments on PEG linkers at room temper- 
ature (2), we employ ksT/Xx = 6pN. Numerically inverting Eq. [3] then 
yields the force extension characteristic fi{li) of a single linker which we 
henceforth will utilize in our model. 

For the pulling geometries studied in our present work, the experimen- 
tally accessible observable is not the magnitude of the total force F = 
y^-_ -i fj, but rather it's normal component with respect to the sample surface 
O, i.e. 

N 

F := F • e, = ^(f, • e,) /, . (4) 

A priori, all quantities in Eqs. [T]|l] are functions of z. However, the actual 
control parameter in the experiment is the displacement 

s = vt (5) 

of the sample surface from its initial position, with a typical experimental 
pulling velocity v = 1000 nm/s. Furthermore, the displacement s can be 



written as the sum of z and the deflection of the cantilever under the force 
F, i.e. 

s = z + FJK , (6) 

with a typical experimental cantilever stiffness (j^ of k = lOpN/nm. Via 
Eq. [6l all the above quantities, in particular the force F from Eq. U may 
thus alternatively be expressed as functions of s, and likewise, via Eq. [5l as 
functions of t. In the following, we will switch between these alternatives 
without much further ado. In particular, F(s), following from Eq. |l]with s 
as independent variable, is called the force extension curve. 

2.5 Rupture probabilities 

For an arbitrary but fixed bond i, the force fi appearing in Eq. H] follows 
by inverting [3] before the bond has ruptured and is given by /» = after 
rupture. According to Evans and Ritchie (jS!), the rupture process itself is a 
thermally activated rate process, governed by 

n,{t) = -k{fi{t))n,{t), mit = 0) = 1 , (7) 

where nj(t) denotes the survival probability of the bond up to time t, and 
k{fi{t)) is the force dependent bond-dissociation rate. Consequently, 



ni{t) = exp 



dt'k{U{t')) 





(8) 



and the probability of bond rupture then follows as —dni[t)/dt. Denoting 
by ti{f) the inverse of the function fi{t), one furthermore obtains the sur- 
vival probability ni{f) := 'ni{ti{f)) of the i-th bond as a function of the 
instantaneous force / acting on that bond. 

In our quantitative calculations below, we will adopt the simplest and 
most common approximation (lid ) 

k{f) = koeMxhf/kBT) , (9) 

where k^ is the force-free dissociation rate and x^ the distance between po- 
tential minimum and barrier along the direction of the applied force. Fur- 
thermore, we will use the following typical values of those rate parameters: 

/i:o = 0.1s"\ Xb/ksT = Q.Iy>^~^ ■ (10) 



We remark that the alternative approximations for k{f) from (jl4| ) did not 
lead to notably different findings than Eq. [9l 



2.6 Simulation and processing of force extension curves 

The initial condition in the form of N bonds at time t = is set according 
to Sect. 12.31 For the further time evolution we adopt - similarly as in 
(2il.l22l.l39l.l4d) - Gillespie's algorithm (45') with time-dependent forces and 



rupture probabilities as detailed in Sects. 12.41 and 12.51 respectively. 

Fig. El exemplifies six force extension curves simulated along these lines. 
They indeed look quite similar to the experimental curves from Figs. [T]A.- 
D. In particular, while in the left panels of Fig. [3] the two rupture events 
("dips" or "jumps") happen to be well separated, the two bonds accidentally 
rupture nearly simultaneously in the right panels. In the middle panels, the 
rupture of the first bond happens to be still visible as a small force dip, but 
due to the limited resolution of a real experimental device (cf. Fig. [ip), it is 
questionable whether in practice such a force signal could be distinguished 
from that of a single bond. 

For comparison, the bottom right panel in Fig. [3] also shows a typical 
force extension curve from a "numerical experiment" with a single bond 
(A^ = 1). This example corresponds to the more likely situation that the 
single bond ruptures earlier than the two bonds, but since these are random 
events with a quite broad distribution (see later sections), it might also be 
the other way round with quite some probability. Within the experimentally 
realistic noise level which has been used for the illustration in Fig. [Sj it thus 
is indeed practically impossible to tell on the basis of the dotted and the 
solid lines in the right and middle panels of Fig. [Sj which one was originally 
due to a double and a single bond, respectively. 

Similar observations apply to more than two bonds. In other words, 
it seems reasonable to assume that the rupture of a multiple bond complex 
cannot be distinguished from that of a single bond if all bonds rupture within 
some small distance Asmax- In view of Figs. [TJ\-D and [3l Asmax = 1 nm 
seems to be an optimistic but still reasonable choice for the resolution limit, 
in goo d agreement with the experimental value of Asmax ~ 4nm reported 
in (jig ) for relatively long linkers. A more detailed discussion of the specific 
choice of Asmax will be provided in Sect. 14.31 

A further consequence of the instrumental noise is that only rupture 
forces beyond some threshold value /min can be clearly distinguished from 
the noise. In our numerical examples below, we will always adopt the ex- 
perimentally realistic choice /min = 20 pN since we found that any other 
choice of /min within reasonable limits (say between pN and 40 pN) hardly 
affected the results with the obvious exception of different "lower cutoffs" 
in rupture force histograms like Figs. [IE^,F. 



As in real experiments, force extension curves which exhibit clearly re- 
solvable signatures of multiple bonds are excluded from the further analysis. 
In turn, this means that force extension curves of multiple bonds will still 
be accepted if all bonds with rupture force larger than /min rupture within 
a distance Asmax- Such multiple bonds, which are erroneously classified as 
single bonds will henceforth be denoted as "false single bonds" or "hidden 
multiple bonds". In turn, both the "true" and "false" single bonds will be 
called "apparent single bonds" . The force value, at which an apparent single 
bond ruptures is defined by the maximum of the force extension curve F{s), 
is termed "rupture force" , and is conventionally denoted by the letter /. Due 
to the random features of tip functionalization (Sect. 12. ip . bond formation 
(Sect. 12. 3p . and bond rupture (Sects. 12.51 and l2.6p . the rupture force / is a 
random variable, whose probability density is henceforth denoted as p{f). 

3 Uniformly distributed receptors 

For a uniform receptor density p^ec (cf- Sect. 12. ip . the number of receptors 
on a surface area ^ is a Poisson distributed random number with mean value 
Prec^- Sampling receptors along these lines, we simulated a number of force 
distance cycles so that 1000 of them were exhibiting at least one rupture 
event (i.e. A^ > 1). In Figs. [1]A,B we moreover averaged over 100 different 
AFM tips (cf. Sect. 12. ip . while Figs. |4lC,D depict representative results for 
one single tip. 

The blue line in Fig. [4]A shows that the probability of classifying an 
actual multiple bond as an apparent single bond is 40-50% and decreases 
only weakly with increasing density prec of receptors. Likewise, the fraction 
of "hidden multiple bonds" among all apparent single bonds (black line in 
Fig. m^) increases remarkably fast as a function of p^ec- However, for these 
Prec values, also the probability to observe at least one rupture event during 
one force distance cycle (red line in Fig. HJ^) is quite large compared to the 
usual experimental findings. We come back to this point in Sec. 14. 1[ 

Fig. SJB shows that the rupture force, averaged over all apparent single 
bonds, exhibits a moderately increasing behavior as a function of the re- 
ceptor density p^ec- This finding is quite plausible in view of the increasing 
number of hidden multiple bonds, which are expected to rupture on the 
average at higher force values than the true single bonds. 

For small receptor densities p^ec, almost only true single bonds contribute 
to the distribution of rupture forces p{f) which is consequently sharply 
peaked and exhibits a fast decay for large / (Fig. HP). For higher receptor 
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densities, the distribution of rupture forces develops a long tail (Fig. [Jp) 
due to the increasing fraction of hidden multiple bonds, much like for the 
real data exemplified by Figs. [IP,F. Interestingly, the distribution still ex- 
hibits a pronounced peak which (by comparison to Fig. [DC) is mainly caused 
by the rupture events of the true single bonds, while no distinct additional 
peaks due to double and more general multiple bonds are visible. 

While Fig. H] was obtained for AFM tips with A^un = 10 linkers. Fig. 
[5] shows that Nn^ = 5 and Nn^ = 15 lead only to minor changes of these 
results. 

4 Clustering of receptors 

4.1 Motivation 

Experimentally, the receptor density is controlled via the chemical prepara- 
tion procedure of the sample surface. Usually it is adjusted such that the 
probability to observe at least one rupture event during one force distance 
cycle (red lines in Figs. H] and [5]) is quite low (typically below 20%). As- 
suming a uniform distribution of receptors on the surface, the probability of 
multiple bonds would then be negligibly small, cf. Figs. H]and[5j It can be 
shown that this observation is largely independent of the remaining model 
parameters and does also apply to other geometries of the AFM tip. Yet, 
in real experiments multiple bonds are actually observed at a much higher 
rate. This seems to us a quite convincing direct evidence that the receptors 
are in fact not uniformly distributed on the surface. Physical reasons of 
why the density of receptors may be locally considerably larger than on the 
average are: (i) The receptors (or the linkers used to immobilize them on 
the substrate) may tend to cluster due to mutual interactions, insufficient 
mixing or substrate inhomogeneities (some areas may be more "reactive" 
than others), (ii) One receptor may exhibit several identical binding sites 
for the ligand, like e.g., streptavidin for biotin (38). 

4.2 Modeling 

Receptor clusters are uniformly distributed on the sample surface according 
to some cluster density pciu- For any given cluster, the number of receptors 
Aciu, is sampled from a Poisson distribution with mean (Aciu) and the recep- 
tors are then independently distributed within a circle of radius i?ciu = 2 nm. 
Other values of -Rciu lead to very similar results, as long as they remain 
smaller than the maximal binding length dmax (cf. Sect. 12. 3p . Likewise, we 
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will always focus on the specific cluster density pc\u = W~^ nm~^ since the 
value of Pciu is in fact irrelevant as long as multiple bonds involving recep- 
tors from different clusters are negligible. Finally, we restrict ourselves to 
average cluster sizes (Ndu) = 2 and discuss the influence of the parameters 
dmax and Asmax in some more detail. 

4.3 Results 

Fig. [6] shows numerical results for different values of the maximal binding 
length dmax- Proceeding as detailed in Sects. [2] and 14.21 we simulated a 
number of force distance cycles so that 10000 of them were exhibiting at 
least one rupture event (i.e. A^ > 1), but unlike in Sect. El no additional 
sampling over different AFM tips was performed. The resulting probabilities 
of observing at least one rupture event during a force distance cycle were 
approximately 5% for dmax = 8nm, 11% for dmax = 12 nm, and 17% for 
'^max = 16 nm. As mentioned in Sect. 14. H these are typical values for real 
experiments. 

Since resolution limits Asmax ^ 1 nm seem to be experimentally unre- 
alistic (cf. Figs. [T]A-D, [3|), Fig. [6K implies that a large fraction of actual 
multiple bonds will be classified as apparent single bonds, and Fig. [6j3 shows 
that a considerable fraction of apparent single bonds are actually multiple 
bonds. 

Also the rupture force distributions p{f) (Fig. [DC) seem to depend only 
quite weakly on the maximal binding length dmax; apart from two details: 
(i) In contrast to Fig. |3JD, the rupture force distributions develop a small 
secondary peak upon decreasing dmax- The reason is that for small dmax, the 
force is always more or less equally distributed among the multiple bonds, 
giving rise to well separated peaks for single, double, triple, ... bonds. With 
increasing (imax values, multiple bonds with unequally distributed forces 
become more likely, hence the force peaks are smeared out and finally dis- 
appear. In either case, it is questionable whether the minima of the rupture 
force distributions p{f) in Fig. [6p can also be resolved in practice, where 
sample sizes are in general much smaller (see also Figs. [IP,F and [7]A-D) 
and where the peaks are broadened by noise, (ii) The complete rupture 
force distribution p{f) for dmax = 16 nm is slightly shifted towards lower 
forces compared to the distribution for dmax = 8 nm. The reason is that the 
rupture force corresponds to the force measured by the AFM cantilever and 
not to the force acting on the bond, see Fig. [2jB. For larger values of dmax 
it is more probable that one is pulling under a large angle relative to the 
surface normal (see Fig. EJB), resulting in smaller values of the measured 



12 



rupture forces. 

4.4 Variations of the pulling velocity 

Fig. [7] shows the numerical imitation of a complete dynamic force spec- 
troscopy experiment: For each pulling velocity v we simulated 500 force dis- 
tance cycles as detailed in the preceding sections, employing dmax = 12 nm, 
^Smax = 1 nm, and the same AFM tip as in Fig. [6j The resulting proba- 
bilities of observing at least one rupture event during a force distance cycle 
were approximately 10% for all pulling velocities v. 

While a detailed quantitative comparison/fitting with any specific ex- 
periment is not the purpose of our present paper, the general similarity 
between the typical experimental data from Figs. [IP,F and the numerical 
simulations in Figs. [7|A-D is quite convincing. In particular, we recover the 
typical "long tails" announced in Sect. [Ij 

For a given pulling velocity v, the survival probability n{f) of the appar- 
ent ("true" or "false") single bonds readily follows from the rupture force 
distribution p{f) according to 



n{f) = j P{f')df' . (11) 

Fig. Ef] presents those survival probabilities for various pulling velocities 
V. As demonstrated in (J2d ). all these functions —v\ap{f) must collapse 
onto a single, ^-independent "master curve" in the absence of false single 
bonds, while experimentally they actually split in a very similar way to the 
curves shown in Fig. [7ji]. In (20 ) this was explained by adopting a para- 
metric randomization of the dissociation rate from Eq. [9l In our case, the 
dissociation rate for a single bond i can be written as A;oexp[(xfeaj)F/A:5T] 
where aj := fi/F. For a fixed force -F, a^ is a random variable depending 
on the number of parallel bonds and their configuration. Effectively, hidden 
multiple bonds thus result in very similar rupture force distributions as bond 
heterogeneities. 

The fact that the curves in Fig. [7f] indeed do collapse onto a single 
master curve for small /-values further corroborates the above discussed 
implications of Figs. |4] C and D, namely that the hidden multiple bonds 
mainly affect the rupture force distribution beyond its maximum, while the 
maximum itself is mainly governed by the "true" single bonds. 
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4.5 Most probable rupture forces 

As seen above, hidden multiple bonds hardly affect the niaxiniuni of the 
rupture force distribution. In other words, the most probable rupture force 
/* admits an adequate and consistent theoretical treatment within the "tra- 
ditional" single bond picture by Evans and Ritchie (fsl), cf. Sect. [TJ In 
particular, the rate parameters Xh/k^T and k^ appearing in Eq. [TU]can still 
be estimated by plotting /* versus In A due to the well established relation 

/* = (kBT/x,) In {x^X/ksTko) (12) 

where A := F'{s{f*))v is the so-called loading rate (jEJ). 

Fig. [7p confirms that the most probable rupture force resulting from 
our simulations indeed depends linearly on In A. By fitting a straight line 
through these points one recovers by means of Eq. [T2]the following estimates 

/i;o = 0.17s"\ Xb/ksT = .0989 p^-^ , (13) 

in good agreement with the original, "true" rate parameters from Eq. [TOj 

5 Discussion 

By a detailed modeling and simulation of a complete dynamic force spec- 
troscopy experiment - including the formation of bonds, their dissociation 
under load, and the post processing of the force extension curves - we have 
shown that multiple bonds cannot be detected with sufficient reliability on 
the basis of the experimentally accessible information, namely force exten- 
sion curves exhibiting several distinct force dips (Figs. [TJ\-D,[3|). In partic- 
ular, in order to explain the typical experimentally observed frequencies of 
force distance cycles exhibiting zero, one, and more than one rupture events, 
we found that assuming some kind of receptor and/or linker clustering seems 
unavoidable (Sect. 14. ip . As a consequence, a quite reliable indicator that 
a significant number of multiple bonds are misinterpreted as single bonds 
is a non-negligible fraction of force extension curves with experimentally 
resolvable multiple dips. 

For the sake of simplicity, we have assumed that the receptors are im- 
mobilized directly on the substrate and that the ligands are connected to 
the AFM tip via linkers of fixed length. Under typical experimental condi- 
tions, about half of the multiple rupture events then turned out to be unde- 
tectable, largely independent of the association dynamics and other details 
of the modeling. Cases when both, receptors and ligands are immobilized 
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via linkers, which may furthermore exhibit a broad length distribution, are 
considerably more difficult to model (|4g) , but intuitively we still expect com- 
parable fractions of hidden multiple bonds, in good agreement with recent 
experimental observations (jlSl ). 

A somewhat related problem in the context of protein unfolding has re- 
cently been addressed by Dietz and Rief (|47l ) . These authors have connected 
the protein of interest to a modular protein chain resulting in characteristic 
sawtooth patterns in the force extension curves. Employing an objective 
pattern recognition algorithm, Dietz and Rief then used these "fingerprints" 
to identify the "true" single molecule unfolding events. Similar experimen- 
tal and theoretical attempts in case of dynamic force spectroscopy are not 
known to the present authors. But at least we can provide a way to "live" 
with those apparently unavoidable "hidden" multiple bonds, namely by fo- 
cusing on the most probable rupture force and disregarding the "long tails" 
of the full rupture force distribution. In other words, our work provides 
a solid justification of what many experimentalists have always been do- 
ing anyway. Finally, Zhu and coworkers (J48) have pointed out already in 
1998 that it is possible to discern single from multiple bonds in micropipette 
measurements just by visually monitoring the detachment of the cell. They 
further showed, that by varying the waiting time for bond formation, it is 
possible to extract the kinetic rates and the mean number of bonds. This 
however required a precise control of the density and distribution of recep- 
tors and ligands, as well as a correct stochastic description of the association 
process. Analogous options in the case of measurements by AFM, where this 
is not available (see Sects. 12.31 and [4.ip . are not known to us. 
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Figure Legends: 

Figure [T} 

(A)-(D) Typical force extension curves from four "single runs" of a dynamic 
force spectroscopy experiment by AFM. Ideally, the force F(s) steadily 
grows with increasing distance s until the chemical bond ruptures (A). In 
(B) and (C) more than one force dip ("downward jump" of F{s)) is clearly 
visible, indicating that more than a single bond was involved. Curves like 
in (D) would commonly still be attributed to a single bond rupture within 
the given noise level and resolution limit. The depicted experimental data 
in (A)-(D) for a protein-DNA complex (PhoB mutant and target sequence, 
pulling velocity of 2000 nm/s, cantilever stiffness 13pN/nm, linker length 
30 nm) have been kindly provided by A. Bieker and D. Anselmetti (Biele- 
feld University). (E) and (F) Typical histograms of experimental rupture 
force distributions for two different pulling velocities v ((TI). After eliminat- 
ing all the experimental force extension curves with clearly visible multiple 
bond signatures as those in (B) and (C), the number of remaining "apparent 
single bond rupture events" contributing to (E) was 202, and 151 for (F). 
The main features are a pronounce first peak (most probable rupture force), 
vague indications of possible secondary peaks, and a "long tail". 

Figure [2j 

(A) Schematic sketch of dynamic force spectroscopy by AFM: a chemical 
bond of interest, e.g. in a ligand-receptor complex, is connected via two 
linker molecules with the tip of an AFM cantilever and a piezoelectric ele- 
ment at distance s. The latter is employed for "pulling down" the attached 
linker molecule at some constant velocity v which in turn leads to an elastic 
reaction force F{s) of the cantilever, determined from the deflection of a 
laser beam. (B) Illustration of the model for multiple parallel bonds. The 
AFM tip is modeled as a half sphere and forces fj act on the ligand-receptor 
bonds. 
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Figure [3} 

The dashed hnes exemphfy six reahzations of force extension curves for 
double bonds (A^ = 2), numerically simulated as described in Sect. [2j 
Receptors were uniformly distributed on the sample surface with density 
Prec = 10~'^nm~^. For a better comparability with the experimental curves 
from Figs. [T]A.-D, we sampled the force extension curves in regular time 
steps of At = 0.1ms and added a Gaussian (thermal) noise with standard 
deviation aj = y/nksT = 6.4 pN (|3). After that, a running average over 
0.5ms was calculated, imitating the effect of an experimental low pass filter, 
and resulting in the solid lines. The distance As between the rupture of the 
two bonds is indicated in each figure. The dotted line in the bottom right 
panel exemplifies one realization of a force extension curve for a single bond 
(A^ = 1, cf. Fig. [DA.)- The somewhat larger fiuctuations observed in Figs. 
[T]A.-D can be attributed to instrumental noise on top of the thermal noise. 

Figure [4) 

(A) Red: Probability of formation (and rupture) of at least one bond (i.e. 
N > 1) within one force distance cycle. Black: Probability of observing 
a "false single bond" (i.e. an apparent single bond is de facto a (hidden) 
multiple bond). Blue: Fraction of "false single bonds" among all multiple 
bonds. All three probabilities are presented for various values of the density 
Prec of uniformly distributed receptors and have been obtained as detailed 
in Sects. [2] and El The error bars indicate the statistical spread (standard 
deviation) due to our sampling of 100 different tips (see main text). (B) 
The corresponding mean rupture forces (/). (C) Representative rupture 
force distribution for one AFM tip (see main text) and prec = 2 • 10~^ nm~^. 
(D) Same for prec = lO^'^nm"^. 



Figure [5j 

Same as Figs. H] (A) and (B), but for Nn^ = 5 linkers in (A) and (B), and 
A^lin = 15 linkers in (C) and (D). 
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Figure [6} 

(A) Fraction of "false single bonds" among all multiple bonds versus res- 
olution limit Asmax for maximum binding lengths dmax = 8nm (solid), 
draa.x = 12 nm (dotted), and dmax = 16 nm (dashed). For further details 
regarding the employed receptor clustering model see main text. (B) Prob- 
ability of observing a "false single bond" (i.e. an apparent single bond 
is de facto a (hidden) multiple bond). (C) Rupture force distributions 
for Asmax = 1 nm. For reasons of better visibility, the distribution for 
c^max = 12 nm is not shown. 



Figure [3 

(A)-(D) Rupture force distributions for different pulling velocities v, assum- 
ing clustering of receptors. For further simulational details see main text. 
(E) The corresponding survival probabilities according to Eq. [TTJ Velocities 
increase in the direction indicated by the arrow. (F) The most probable 
rupture force /* from (A)-(D) versus logarithm of the loading rate A. The 
solid line shows the best linear fit. For more details see main text. 
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